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Gravitational waves in an anomaly-induced inflation 
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The behaviour of gravitational waves in the anomaly-induced inflationary phase is studied. The metric per- 
turbations exhibit a stable behaviour, with a very moderate growth in the amplitude of the waves. The spectral 
indice is computed, revealing an almost flat spectrum. 



The goal of the present study is to analyze 
the fate of gravitational waves in a background 
defined by an inflationary scenario created by a 
trace anomaly induced by quantum effects gen- 
erated by matter fields in the primordial Uni- 
verse [(T). This model is a generalization of the 
Starobinsky model [ 2^, where a similar mecha- 
nism has been developed. The counterterms nec- 
essary to avoid divergences in the quantization of 
these conformal quantum fields in a Friedmann- 
Robertson- Walker background leads to a higher 
derivative action with a trace anomaly. Details 
of the model may be found in [ 1 . An important 
point concerns the fact that an inflationary phase 
may be obtained with this effective gravitational 
action. When all fields are massless, this infla- 
tionary phase is an eternal de Sitter phase, where 
many of the traditional problems of the inflation- 
ary scenario appear. However, if massive confor- 
mal fields arc included, the inflationary phase has 
more complicated form. A detailed analysis of the 
massive case reveals that a huge amplification of 
the scale factor may be obtained, of the order 
of 10^° e-folds. However, from the observational 
point of view, only the last 65 — 70 e-folds are rel- 
evant. During this final period, the behaviour of 
the scale factor may be approximated by an ex- 
ponential function, the Hubble parameter being 
constant. It is possible that a transition to the 
FRW standard scenario may be achieved solving 
the graceful exit problem [i4j. 

The equation governing the behaviour of grav- 
itational waves in the anomaly-induced model, in 



terms of the cosmic time is 
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bo = bo, 6i = Hbi , 62 - H% , 
bs = H% , 64 = H% , 

ni ^ Hhi , = ^2 , «3 = "3 ■ (2) 
Here the quantities with tildes are pure numbers: 



60 = 1, ^1 = 6, 62 = 11, 63 = 6, 
rii = -2 , n2 = -2 , = 1 . 



(3) 



while 64 depends on the multiplet content of the 
matter fields: 
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where A'o, A^i/2 and A^i are the scalar, fermionic 
and vectorial numbers of matter fields. For ex- 
ample, in the minimal standard model, A^o ~ 4, 
A^i/2 = 24 and A^i = 12, leading to 64 ~ 0.2. 

As usual, H = a/a is the Hubble parameter. 
The parameters described above take those value 
in the last 65 e-fold. In this case, the background 
model enters in a quasi-de Sitter and the scale 
factor behaves essentially as 
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where, in an appropriate GUT model, 
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Changing to the conformal time 77 such that 
dt = adrj, the scale factor takes the form 

1 



(6) 



with ?7 < 0. The Universe expands as ry 0_. 
In terms of this new time parameter, and per- 
forming a Fourier decomposition of the function 
h — h{x,rj) — h{T])e^'''^, k being the wavenum- 
ber of the perturbation, we obtain the following 
fourth order differential equation describing the 
evolution of gravitational waves. 



/,» + 2fc2/j"+|^ + fc4|^^0. (7) 

The units are such that k = I implies a comov- 
ing wavelength of the perturbation of the order 
of the Planck's length. Let us notice the remark- 
ably simple form of the equation iQ), which is 
essentially simpler that the equivalent equation 
in terms of the physical time ^ . 

We first investigate the behaviour of h in 
the two asymptotic regimes: small and large 
wavenumbers. For large wavenumbers A: >> 1 
and we may approximate the equation to 

/i*" + 2k^h" -I- fc4/i = , (8) 

whose solution is 

h = c±e±''^''' + c'^r^e^'^-"! . (9) 

Since 77 approaches zero as the universe expands, 
the solution above is a combination of oscilla- 
tory and decreasing oscillatory modes. A quan- 
tum state can be implemented in the initial state 
(when 77 — !■ —00), if, for example, c_ = co/V2k, 
Co = constant, c+ = c'j_ = 0. For small wavenum- 
bers, on the other hand, the equation simplifies 
to 

/i» + %/i = 0. (10) 
T 

We can look for solutions under the form h oc rjP, 
where p is a number. In this case, p obeys the 
algebraic relation 



There is no root for p with Rep < 0. Hence, there 
is only decreasing modes in the long wavelength 
approximation. 

All the behaviours sketched above were con- 
firmed through numerical integration. The re- 
sults indicate that the inflationary phase in this 
model is quite stable with respect to tensorial per- 
turbations. There is a small production of gravi- 
tons during the inflationary phase. This seems to 
be a positive result since a very large production 
of graviton during the inflationary phase would 
lead to undesirable consequences since a large am- 
plification of gravitational waves brings the prob- 
lem of back reaction [ 5^ and, at same time, ren- 
ders the linear approximation doubtful. No such 
problems exist in the behaviour of gravitational 
waves in the context of the anomaly-induced in- 
flationary scenario. 

It is instructive to compare the model described 
above with the results obtained with the case of 
the standard inflation based on a cosmological 
constant term (or an evolving scalar field) [ ^ . 
Since a de Sitter phase is a good approximation 
to the standard case, the scale factor behaves as 
before. In this case the equation governing the 
behaviour of gravitational waves takes the form 
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The solution for this equation is 



h=^'i ciH'^y^ikr,) + C2H^f^{k7]) 



h = 0. (12) 



(13) 
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where the Hi, ' are Hankel functions of first and 
second kind, and the c's are constants. 

From the point of view of gravity production, 
directly connected with the amplification of the 
gravitational waves, the anomaly-induced infla- 
tionary model is very different from the standard 
one. While the solutions H13|l exhibit an ampli- 
fication of the order 10^^ in the last 65 e-folds, 
in the anomaly-induced infiation the amplitude of 
the gravitational waves remains essentially stable. 
The huge amplification of gravitational waves in 
the standard inflationary scenario brings essen- 
tially two problems: first, the non- linear regime 
is reached very quickly unless the initial pertur- 



p(p-l)(p-2)(p-3) + 64 = 0. 
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bations are excessively small; moreover, it is difS- 
cult to ignore in this situation the back-reaction. 
The firsts problem may be coped with through 
the quantum normalization: the initial ampli- 
tude must be of the order of the inverse Planck's 
mass. The second one is object of many inves- 
tigations today [ 5 . These problems are absent 
in the anomaly-induced inflationary model: since 
a perturbation originated in the beginning of the 
de Sitter phase would have been amplified at the 
end of this phase by a factor of order of unity, the 
perturbation enters in the radiative phase very 
small. During the radiative phase, the equation 
governing the evolution of gravitational waves is 
given by H12|l . with a oa i] and with a cx 7]^ during 
the matter dominated phase. The gravitational 
waves will be considerably amplified only during 
these phases. The non linear regime may be at- 
tained only quite recently. 

We turn now to the problem of evaluation 
of the spectral indice of the power spectrum in 
the anomaly-induced inflationary scenario. The 
power spectrum is defined as 

P'{k)^k'Sl. (14) 

The spectral indice is obtained through the rela- 
tion 
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where the subscript T indicates that we are evalu- 
ating the power spectrum for the tensorial modes. 
A scale invariant spectrum is characterized by the 
fact that the amplitude of the perturbation are 
the same for any value of k when the perturba- 
tion cross the horizon. This implies nx — 0, when 
the expression H14|l is evaluated at the horizon 
crossing. 

A fundamental point in the analysis to be per- 
formed is the question of initial conditions. In 
the anomaly-induced inflation, as in the standard 
inflaton-based inflation, the seeds of the initial 
perturbations are quantum mechanical. The nor- 
malized initial quantum spectrum takes the form 

p-ikr] 

Mry),^^ . (16) 



2k 

Let us consider the initial time rji ^ 1 in 
Planck's unities. The final time is 77/ ^ 10^"^°. 



Suppose also that the inflationary phase ends at 
about t ~ 10"'^* s, consistent with the fact that 
Hq ^ lO^^Mpj. Since in our unities fc '--^ 1 
means a initial perturbation of the order of the 
Planck's scale, perturbations today in range from 
the Hubble horizon and some hundreds of mega- 
parsecs, evaluated today, implies initial pertur- 
bations with 27rlO~^ < fc < 27rl0"^ We now 
compute numerically the perturbations using the 
initial conditions (|16|l for the traditional inflation- 
ary scenario and for the anomaly-induced infla- 
tion. In both cases, the spectrum is essentially 
flat: ht ~ at the end of the inflationary phase. 
Considering, for example, tt/IOO < k < tt/IO, 
the spectral indice is ut ~ —0.027. This may be 
compared with the spectral indice for the stan- 
dard inflationary scenario which, for a de Sitter 
phase, is strictly zero in the large wavenumber 
approximation. 
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